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Abstract 



A new class of lattice Dirac operators D which satisfy the index theorem have 
been recently proposed on the basis of the algebraic relation 75(75!)) + {'y^D)^^ = 
Q^i 2a^^^^(75-D)^'^''~^. Here k stands for a non-negative integer and k = corresponds 

I to the ordinary Ginsparg- Wilson relation. We analyze the locality properties of 

■ Dirac operators which solve the above algebraic relation. We first show that the 
free fermion operator is analytic in the entire Brillouin zone for a suitable choice 
of parameters mo and r, and there exists a well-defined "mass gap" in momentum 

^ ' space, which in turn leads to the exponential decay of the operator in coordinate 

space for any finite k. This mass gap in the free fermion operator suggests that the 
operator is local for sufficiently weak background gauge fields. We in fact establish 
fSJ I a finite locality domain of gauge field strength for r5 = 75 — (a^^D)'^^'^^ for any 

■ finite k, which is sufficient for the cohomological analyses of chiral gauge theory. We 
also present a crude estimate of the localization length defined by an exponential 
decay of the Dirac operator, which turns out to be much shorter than the one given 

^ ' by the general Legendre expansion. 

&■ 

^ ■ 1 Introduction 



We have recently witnessed a remarkable progress in the treatment of lattice fermions[l]- 



' [4]. The first breakthrough may be traced to the domain- wall fermion [5] [6], which was 

^ followed by the overlap fermion[7]. See also related works in [8]. It is well known that 

the overlap fermion[2], which developed independently of the Ginsparg- Wilson relation[l], 
satisfies the simplest version of the Ginsparg- Wilson relation. The recognition that the 
fermion operator, which satisfies the simplest version of the Ginsparg- Wilson relation, 
gives rise to the index theorem on the lattice [3] and thus modified but exact lattice 
chiral symmetry [4], was crucial in the recent developments. The locality properties of the 
Neuberger's overlap operator have also been established [9] [10]; the operator is not ultra- 
local [11] but exponentially local, which is considered to be sufficient to ensure locality in 
the continuum limit with the lattice spacing a — > 0. 

In the mean time, a new class of lattice Dirac operators D have been proposed on the 
basis of the algebraic relation[12] 

75(75/^) + (75^5)75 = 2a''+\j,DY''+' (1.1) 



1 



where k stands for a non-negative integer, and /c = corresponds to the ordinary Ginsparg- 
Wilson relation[4] for which an exphcit example of the operator free of species doubhng 
is known [2]. It has been shown in [12] that we can construct the lattice Dirac operator, 
which is free of species doublers and satisfies a lattice version of index theorem [3] [4], for 
any finite k. An explicit calculation of the chiral anomaly for these operators has also 
been performed [13]. Here 75 is a hermitian chiral Dirac matrix and 75D is hermitian. 
In a sense, this new class of fermion operators are regarded as the first generation of 
lattice fermion operators which are directly motivated by the Ginsparg- Wilson relation. 
In contrast to the most general form of the Ginsparg- Wilson relation [1], our algebra 
(1.1) allows an explicit solution which satisfies the index theorem: Our explicit solutions 
illustrate the features of possible solutions of the general Ginsparg- Wilson relation. 

A salient feature of Dirac operators corresponding to larger k is that the chiral sym- 
metry breaking term becomes more irrelevant in the sense of Wilsonian renormalization 
group. For H = a'j^D in the near continuum configurations we have, for example. 



respectively. The first terms in these expressions stand for the leading terms in chiral 
symmetric terms, and the second terms in these expressions stand for the leading terms 
in chiral symmetry breaking terms, respectively. This shows that one can improve the 
chiral symmetry for larger k, though the operator spreads over more lattice points for 
larger k. As another manifestation of this property, the spectrum of the operators with 
/c > is closer to that of the continuum operator in the sense that the small eigenvalues 
of D accumulate along the imaginary axis ( which is a result of taking a 2A; -|- 1-th root), 
compared to the standard overlap operator for which the eigenvalues of D draw a perfect 
circle in the complex eigenvalue plane. 

In this paper, we analyze the locality properties of this general class of lattice Dirac 
opeartors. As for the overlap Dirac operator[2], which corresponds to /c = 0, the very 
detailed analyses of locality properties have been performed by Hernandez, Jansen and 
Liischer[9], and Neuberger[10]. A similar locality analysis of the domain wall fermion has 
been given by Kikukawa [14]. We establish the locality properties of the general class of 
operators (1.1) for all finite k, which are sufficient for the cohomological analyses of chiral 
gauge theory [15] [16], for example. 

2 A Brief Summary of the Model and Notation 

The exphcit construction of the operator, which satisfies the relation (1.1), proceeds by 
first defining 



H ~ 



H ~ 



75ai p + 75(750? P)^ for A; = 0, 
75ai p + 75(750? P)^ for k^l 



(1.2) 




(2.1) 
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The operator D^*^"*"^^ is in turn expressed as a generalization of the ordinary Wilson Dirac 
operator as 

Dg'^+i) = + (By>^+' - (H^y^+K (2.2) 

a 

The ordinary Wilson Dirac operator D^, which corresponds to D^) , is given by 

Dwix, y) = i'j'^Cf^ix, y) + B{x, y) - -mo5^^y, 

a 

T 

za ^ 

Uf,{y) = exp[iagA^{y)], (2.3) 

where we added a constant mass term to Dw- Our matrix convention is that are 
anti-hermitian, (7^^)^ = —7^, and thus = '-)^C^{n,m) is hermitian 

(2.4) 

To avoid the appearance of species doublers, we need to satisfy 

2r > mo > (2.5) 

and the choice 2rn^^^ — 1 simplifies the normalization of the operaror D in the continuum 
limit as in (1.2). 

By defining hermitian H = a'j^D, the algebra (1.1) is written as 

75// + i/75 = 2i/2fe+2 (2.6) 
or equivalently as a set of equations 

H^^, - ^,H^ = 0. (2.7) 

The first of these relations, when regarded as an equation for is solved by i?(2fe+i) 
in (2.1). 

The solution of the original algebra (2.6) is then given by 

H = (i/(2,+i))V(2'=+^) (2.8) 

in the representation where H(2k+i) is diagonal and the value of H is chosen as real. This 
construction means det H — (det Hi^2k+i)y^^'^''~^^^ ■ The case A; = is reduced to the overlap 
Dirac operator. 

When one defines 

Ts = 75 - H^2k+i) (2.9) 

the defining algebra is written as 

r^H + HF^^O. (2.10) 
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The (fermionic) index is then given by 



Trr5 = n+-n_ (2.11) 

where n± stand for the the number of zero modes of 

H^n = (2.12) 

with 75(^n = ^fn, respectively. In the smooth continuum hmit, one can estabhsh the 
index theorem 

2 

n+-n^^l ci^x J^tre'^'^^^F^.F,^ (2.13) 

for all (finite) k. The right-hand side of this relation stands for the Pontryagin (or instan- 
ton) number. See Ref.[12] for further details. 

To analyze the locality, one thus has to examine both of 

H^2k+i) ^ H"'+\ (2.14) 

which appears in Fs, and H itself. For the free fermion case, where one can write the 
operator H explicitly, it turns out that the localization length is essentially the same for 
both of H and i?(2fc+i)- We later argue that this property holds for operators H and 
H{2k+i) with non-trivial gauge fields also. 

3 Mass Gap for Free Fermion H^^^^ 

For vanishing gauge fields, we have 

Y^D^{x,y)e^^y = [-^7'^^^ + ^^(l-cosafc,)-^K'=^ 



Dw{ak^)e''^ (3.1) 



and thus 



where we defined 

= afcp,. (3.3) 

Note that we are working in the infinite lattice space, and thus we have a continuous 
spectrum in the Brillouin zone, which is chosen to be 

TT 37r , , 

-^<P.<^ (3.4) 
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for all /X = 1 ~ 4. 

Similarly we can define D^r''~^^\p), or 



-{mof'^'. (3.5) 

We then have 

The general operator in the coordinate representation is given by 

i/(2.+i)(x, y) = ^ / ^^(2fe+i)(PM)e^^(^-^)/". (3.7) 
In the denominator of the general class of Dirac operators H{2k+i) (Pn) , we thus have 

= i^T^' + {[E(i - c.)]''^' - (^0)^^=+^}^ 

= [E(i-4)]'^-^^ + {[E(i-CM)]''=^'-K)^'=+^f (3.8) 

where we defined the variables 

= sinp^, 

= cosp^ (3.9) 
and we use the value of the Wilson parameter 

r = l (3.10) 
throughout the present paper. We also set 

ml''+^ = 1. (3.11) 
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3.1 Minimum of 

To analyze the "mass gap" of we examine 

= [E(i - cDf'^' + eV - c,)T^' - 2E(i - c,)r^' 

+[E(1 - c.)r^''^'^ - E(l - c,)T^' + 1. (3.12) 

We first observe that 

E(i - s)]'^''^^^^ - E(i - > (3.13) 

by noting 

E(l - Cm)]' > E(l - c.)'] > (3.14) 

where the equality holds only for = 1 for d — 1 indices /i in d>2. 
We next define 

^ = E(1-Cm)>0, 

A + S = E(l-c^)^ (3.15) 

then one can confirm 

->l + S = -(A + S) + 2S = E(l-cJ). (3.16) 



We thus have 



E(i - cDr^' + E(i - c,)T^' - 2E(i - c,)] 

n n n 

^ (_A + Bf''+^ + {A + Bf^-^^ - 23'"'+^ 

I 



2k+l 



2 E (^^t^]A'B^'+'-'>0 (3.17) 



l=even,^0 ^ 

since B > 0. Here the equality sign holds only for 5 = A = with = 1 for all or 
A = but S ^ with C/,(l - c^) = 0. 

By this way, we have established that for any finite k 

F(k) > 1 (3.18) 

where the equality holds only for = 1 for all n or one of = and the rests of = 1 
in the space-time dimensions d > 2. Namely 

||^f+^)||^||75^£'^+^^(p,)||>l (3.19) 
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for niQ = 1. If one combines this relation with 

^ = s|Tr(^.^^ir'V) = W,75^) (3.20) 

where A(mo) is the eigenvalue of the hermitian operator H^!'~^^\ we obtain a "flow in- 
equality" [10] 

l^^l = l(</',75</')|< 110111175011 = 1. (3.21) 

We can then establish the "mass gap" of ||i^|^^^^''(mo)|| for other choices of ml^~^^. Fol- 
lowing the analysis in Ref.[10], one can establish, for example^], 

\\m'^''\\^\h.m'^'\p,)\\>l (3.22) 

for 

2ml^+^ = 1 (3.23) 
which simplifies the normalization of the continuum limit as in (1.2). 

3.2 Maximum of 

To analyze the maximum of we define 

F(,)(c) = [l-c' + + {[1 - c + S^if^+i - ly (3.24) 
where we defined c as one of four variables c^, and 

s(^) = E(i - 4) > 0, 

S« = ^(1 - c,) > 0. (3.25) 



Here, the summation runs over the three indices /i except for the one used for c = c^. We 
then obtain 

F(',)(c) = -2(2A; + 1){[1 - + S(2)]2*^c + {[1 -c+E^'^^' - -c+E^'^}. (3.26) 
We now examine one special section of species doublers specified by 

c^^<0 for all fi. (3.27) 



^In this case, the equahty holds only for = 1 for all fi. 
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In this domain, the above formula of -F(j.)(c) gives rise to 

= {[1 - + S('^]''c + {[1 - c + _ i}[i _ c + ^'j 

> [1 - c + E(if ^={0 + [1 - c + S(if ^=+1 - 1} 

> [l-c + E(if'={[l-c + EW]2'=+i-2} 

> [1 -c + E(i)]2'={[4]2'^+^ -2} > 0. (3.28) 

Here we used 

0<i-c2 + E(2) = (i_c)(l + c) + ^(l-c^)(l + c^) 

< 1-c + eW (3.29) 

in the above domain, and thus 

[1 - + yS% > [1 - c + E«]c. (3.30) 

This shows that F'q,-^{c) < for any choice of c = c^. Thus the maximum of F(^k) appears 
at = — 1 for all /i in this domain. 

This value is also the absolute maximum of F(^k) in the entire Brillouin zone. This is 
because 

F(.)(ci,c2,c3,c4) = [E(i-^)r^'+{[E(i-CMF'^'-ir 

> -^(fe)(-Cl,C2,C3,C4) (3.31) 

if all Cn < 0: One can in fact confirm that F(fc) (ci, C2, C3, C4) is not smaller than the function 
obtained from it by reversing the signatures of some of c^. This means that the maximum 
of in the domain < for all // is in fact the absolute maximum of F(^k) -i which is 
given by setting all Cn — —1 

Fik) < (8^*^+^ - 1)^ (3.32) 

The absolute maximum of F(fe) = {D^!'~^^\pf^)y D^!'^^\p^) provides the upper bound for 
the eigenvalue and thus the norm of the operator (D^J'^' V^!J'^'^ ■ 

3.3 Locality of Free H(2k+i)ix,y) 

Following the analysis of Ref.[9], we can then establish the locality bound for the free 
operator D(^2k+i){x,y) = 75-H'(2fe+i) (2;, y) in (2.1) by 

1 K , 9\\x — vW , , , 

<--exp[- , . ' ]■ (3.33) 



(^i^'+'))tDi^'=+') 1-^ 2(2A; + l)a^ 
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When one denotes the maximum and minimum of by v^ax a-nd Vmin, respectively, the 
parameters are defined by 



cosh 6* 



t = e' 



For Vmax > Vmin, wc obtain 



(3.34) 



(3.35) 



and the localization length, which is defined by the exponential decay ~ exp [— | |a; — |/ 1 1 /L] , 
is estimated by 



In the present case 



One thus obtains 



or 



Lc^{2k+l)aJ^. 

(82^+1-1)2 + 1 

cosh B = ' , , , -; 

(82fc+l _ 1)2 _ 1 

cosh^ ~ 1 + 2 X 8-2{2fe+i) 



^ ~ 2 X 8-(2fe+i). 
The localization length L is then estimated by 

L~(2/t + l)23(2^+^)a 

or in 2-dimensional case 

L ~ (2A; + l)22(2fc+i)a. 



(3.36) 

(3.37) 

(3.38) 
(3.39) 

(3.40) 
(3.41) 



We will later show that this localization length is much bigger than our estimate on the 
basis of analyticity, which is closer to the numerical estimate. 

We here present a suggestive interpretation of the above estimate of cosh^ (3.34). The 
above formula is written as 



coshe = l+ 



This last relation is re-written as 



COS TT — COS 



;cos(0 + i^) -cosO) = 



(3.42) 



(3.43) 



by noting cosi^ = cosh^. This formula is interpreted that one is estimating the zero of a 
generic function f{cos9) of a single variable cos^, which has a minimum f{cos9) — Vmin 
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at cosO = 1 and a maximum f{cos9) — Vmax a-t costt — —1. As a crude approximation 
for such a function, one may consider a function 

/(cos^) = Vrnin + (cos g - cosO) (3.44) 

COS TT — COS 

which has an average slope 

COS TT — cos 

at cosO = 1. One may then look for the position of the zero of /(cos 6*) = 0, which is 
located closest to the real axis. This equation corresponds to (3.43) above. The zero of 



/(cos 6') gives rise to a singularity in y'/(cos6'), which in turn gives rise to the localization 
length in an analogous analysis in Section 5 later. 

This estimate of the zero (3.43) may be a reasonable one if one knows only the values 
of Vmax and Vmin in the entire domain — 1 < cos^^ < 1 but no more information about 
the functional form of f{cos6): But this estimate becomes poor if all the derivatives of 
/(cos 6*) at cos 6* = 1 vanish f^'-^cosO = 1) = up to a certain integer /, for example. 
More about this interpretation will be commented on later. 

4 Free Fermion Operator H = 075!) 

We first define an explicit form of H in the free fermion case. We start with an ansatz of 
the general solution as (see also Ref.[17]) 

H^^^,[A + B^] (4.1) 



and evaluate 2k + 1-th power. We first note 

= \[A - B^][A + B^] = \[A'' + 5^4] (4-2) 
which satisfies the necessary condition (2.7) 

75//' - //'75 = 0. (4.3) 

Here, we defined 

= ^(sinaA;^)2 > 0. (4.4) 

Our convention is (7^^)^ = —7'*. Our ansatz thus gives 

jj2k+i ^ (l)2fc+i[^2 ^ B^^]''j5[A + B'^] (4.5) 
2 d Qj 
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On the other hand, our exphcit construction uses the building block 

a^^ a ^ a a 



and 



We now define 

A*' 

(X 

^ E-(l-cosa^,r+i-(^r^+\ 
M ^ a 

Then our construction gives 



Comparing these two expressions (4.5) and (4.9), we have the relations 



2 0? 



From this we obtain 

by noting the relation 
and thus 



2 « 



LL 



By this way we obtain 



A = 22fe+l 1 H ; 2fc+l 
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k 1 A;+l j fc + 1 



k 1 ,S^,I,r Ml. , fc 



2 



fe ^ 1 ^ fc+i 



fe 



Namely 



-22W+i^^=)w+i^^H^-Mk)—K (4.14) 



1 fc+1 1 fc+1 / fc+l / fc ^ 

H = 75(:^)^(^=)^{(V^^ + Mfc)^-(7i/2,-Mfe)^^} 



, 1 , fc+1 , 1 , 



= 75(^)2'=+i(^=)^'=+i{(V^(fe) + ^fc)2'=+i - (y^F(fe) -Mfc)2fe+i ^ 



where 



(4.15) 



Mk = [E(l-CM)r+'-l- (4-16) 



We have estabhshed that > 1 in Section 3. To estabhsh the analyticity of H in 
the entire Brillouin zone, we thus have to examine 



f± = /^±Mfe 



= Jisr'^' + {[E(i - c,)]"^^' - 1}' ± {[E(i - c.)]''^' - !}• 

(4.17) 

If we find a non-analytic behavior of 

or (4.18) 

inside the Brillouin zone, we are in danger to encounter the non-locality of our operators 
H in (4.15). 



We first note that we can estabhsh 



f->l (4.19) 
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in the physical domain by using 

/- = iE(i - ciw'^' + {[E(i - c,)?'^' - - [E(i - 

+ 1 

[E,(i-<)]^^+^-[E,(i-c,r+^+i 



1 

+2 

> i (4.20) 

since [Z^^(l — c^)]^*^"^^ — [E/^(l ~ c^)]^*^''"''^ > in the physical domain. We also used the 
gap relation ^[E;.(l - cl)]^'+' + {[^^{1 - c,)]^'+' - 1}' > 1- 



We next note the identities with a non-negative integer k 



/f +^ = (s^)'=+V/i'=+^ = - cDf+'/fl"^' . (4.21) 

fc 

This shows that we can establish the analyticity of in the entire Brillouin zone, if 

we can show that /+ is positive and non- vanishing in the domains of "species doublers" . 
In a domain of species doublers, for example, ci < and other > 0, we have 



\ 



[1 - + E(i - ^W^' + {[1 - + E(i - - 1}' 

+[i-ci+E(i-CM)r^'-i 



> 



[1 - E(i - + {[1 - + E(i - - 1}' 

> 1 (4.22) 

where the equality does not hold simultaneously. Here Yl'^ stands for the sum over = 
2 ~ 4. We thus have 

/+ > 1. (4.23) 

By repeating a similar analysis of /+ for other domains of the species doublers, we can 
establish the positivity of /+ in all the domains of species doublers and thus the analyticity 

k fc+1 

of /^''+^ in the entire Brillouin zone. Similarly, we can establish the analyticity of f^^^ 
in the entire Brillouin zone, by using the above identity (4.21) and the positivity of /_ in 
the physical domain. 

Wc have thus established the analyticity of the free H for any finite k in the entire 
Brillouin zone. 
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5 Locality with Exponential Decay 



Wc analyze the asymptotic behavior of the operator D by means of the Fourier transform 

k 

of the momentum representation. We first examine /_ (c^) ^fc+r appearing in H with 



-[E(i - c.)]''^' + 1 (5-1) 

as a function of complex by fixing the phase convention in the physical domain. It can 
be confirmed that /_ vanishes only for = 1 for all even if one analytically extends 

k 

f- as a function of c^^. But at c^^ = 1 for all /x, /_ = 2 7^ and thus /_(c^)2fe+T is regular 
in the neighborhood of the domain of physical species. One can also establish that /+ is 
non- vanishing in the neighborhood of the domain of species doublers even if one extends 

k 

analytically in c^. By recalling the identities (4.21), one can define /_ (0^)2^+1 in the 

k 

domain of species doublers by using /+. By this way, one can establish that f-{c^) 2^+1 is 

k 

analytic in the neighborhood of the entire Brillouin zone. To be more precise, /_ (c^u) 2^+1 
is analytic in the strip of complex plane with = cosp^, 

—00 < '^p^ < 00 (5.2) 
except for the singularities appearing in 



Fk{c,) = /(E(l - clW'^' + {E(l - c,)?'^' - (5-3) 



Similarly, one can show the analyticity of /+(c^)2'!+i in the same domain of complex p^ 



except for the singularities of y'Ffc(c^). One thus establishes that the numerator factors 
of the free fermion operator H in (4.15) do not induce singularities other than those of 
Ffc(c^) even in the complex plane. 

We thus encounter the singularities of the operator H only at the singularities of 
Fk{c^j)-i which appears in the denominator of in the complex plane. The existence 
of the mass gap means that those singularities are away from the real axis by a finite 
distance. It is important to recognize that the possible singularities are located at the 
same points for both of free H and free H'^^'^^ — H{2k+i)i which are controlled by the 



behavior of yFk{cn). Since we have shown that free H^''+^ — H(2k+i) is local for any 
finite k in Section 4, we infer the locality of free H for any finite k. 



In the following we establish the locality of free H directly by examining 

1 /■37r/2 
t/2 



D{x^, 0) = - l^H{p^) eMw,.x^/a]d^pJa^ (5.4) 
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Impi 



X 



n 



X 



7r/2 



37r/2 Repi 



Figure 1: We can shift the integration path of pi upto the singularity closest to the real 
axis. Crosses in the figure represent the singularities. The line integrals along the vertical 
lines at ^pi = 37r/2 and ^pi — — 7r/2 cancel. 



for large positive x^. Our metric convention is g^j^^, = (1, 1, 1, 1). We replace 



in the numerator of H{p^), and thus this factor does not influence the analysis of locality. 

We set x'' = for /i 7^ 1. 

When one lets positive and very large, we shift the integration path of pi to a line 
parallel to the real axis in the upper half complex plane. We note that the integrand 
75if(p) is identical at pi = —tt/2 + iy and pi = 3n/2 + iy for a real y. Also, x^ assumes 
the values only on the lattice points and thus 



The contributions in the line integral along the vertical lines at 3?pi — 3Tr/2 and 3?pi = 
— 7r/2 thus cancel. 

The leading behavior in x^ is thus determined by the integration path which is parallel 
to the real axis and flrst touches the singularity (i.e., the singularity closest to the real 
axis) in the complex pi plane(See Figure 1). See also Ref.[18]. 

This property is also confirmed directly as follows: We flrst recall the Riemann- 
Lebesgue lemma for a regular H{pij) 



and consider the case where the singularity closest to the real axis is written with real gi 
and positive ri as 



■/^ ia Q 



(5.5) 



exp[i(37r/2 + %y)x^ lo\ = exp[i(— 7r/2 + iy)x^ /a]. 



(5.6) 



lim D(x'^, 0) = lim 




Pi = qi + iri. 



(5.8) 
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We then consider the integration path parallel to the real axis given by 



Pi = Qi + i{ri - e) (5.9) 
where e is a small positive number. The above integral is then written as 

D{x^,0) = exp[-(ri -e)a;Va] (5.10) 

1 /■37r/2 ' 

X - / l^H{p^) exp[igixVa + Y. ip^x^" / a]dqid^p^/ a^. 
a J~n/2 y 

Since the integrand is still regular, one can apply the Riemann-Lebesgue lemma to this 
last integral. This means the exponential decay of the operator with (a conservative 
estimate of) the locahzation length L given by 

L = aj{ri-e). (5.11) 

The existence of the well-defined mass gap thus leads to the exponential decay of the 
operator for a large spatial separation. If the singularity of H{p) in pi, which is closest 
to the real axis, is written as pi = gi + iri, the above analysis gives an estimate of the 
localization length L at 

L^a/n. (5.12) 

Moreover, this estimate of the localization length is the same for H and for if^fe+i — 
H{2k+i), since we have shown that the singularities in the complex pi (or ci) plane are 
identical for these two operators. 

In the one-dimensional integral, one has 

/oo 
dyexp[—\x — y\/L] exp[—\y — z\/L] = (L + \x — z\) exp[— |x — z\/L] (5.13) 
-oo 

This shows that a multiplication of two operators, which decay exponentially, produces 
an operator which decays with the same exponential factor up to a polynomial prcfactor. 
A generalization of this relation suggests that a suitable 2/c + 1-th root of an exponentially 
decaying operator gives rise to an operator with an identical localization length for any 
finite k. This is what our analysis of the free fermion operator so far indicates. Later, we 
establish a finite locality domain of background gauge fields for the operator i^^fc+i for 
any finite k. Although we cannot explicitly demonstrate the locality domain of H in the 
interacting case, we infer the locality domain of H from the above analysis of free fermion 
operators and a generalization of (5.13). 

Incidentally, if the operator H{p^) is analytic in the entire strip in the complex plane, 
i.e., if ri = oo in (5.8), the localization length vanishes and the operator becomes an 
ultra-local one. 

5.1 Crude Estimate of Localization Parameters 

We here present an estimate of the position of singularities in H, which is closest to 
the real axis, and illustrate an estimate of the localization length. The purposes of the 
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analysis of a simplified version in this subsection are two-fold; firstly, we illustrate the 
explicit evaluation of the localization length, and secondly, we show that the localization 
length can generally be much shorter than the estimate given by the general Legendre 
expansion in [9]. We examine the position of zeroes in 

= [i_ci]2(2fc+i) + i + [i_c2]2'=+i-2[l-ci]2'=+i = (5.14) 

This corresponds to the simplest d = 1 dimensional example, which is obtained from the 
general formula by setting = 1 (/i 7^ 1). Since the values of the variables = 1 
(yU 7^ 1) in (5.14) give the stationary point of the integrand of the Fourier transform of 
our problem with = (/i 7^ 1), this example may also have some implications on the 
case d > 1 for k ^ 1: The stationary points have a non- vanishing phase space for the 
Fourier integral and thus they are expected to give an important contribution. 

The singularity in ci closest to the real axis gives a dominant contribution for large 
Xi. We thus look for the solution of 

f,{z) = + 1 + z2fe+i[(2 - ^)2fc+i - 2] = (5.15) 

where we set z = 1 — ci . 

The simplest case k = gives rise to0 

fo{z) = 1 (5.16) 

and thus no singularity with ri = cxd, namely an ultra-local operator. This estimate agrees 
with the one on the basis of the Legendre expansion 

and thus 

^^^j^^^^W^tW (5.18) 
For A; ^ 1, we assume that ^^fc+i jg g^all for the solution. Then 

, ^ (-1)21^1 = exp[±^^]i (5.19) 

where we chose the solutions closest to the real axis in terms of pi defined by Ci = cospi. 
This value gives 

(5.20) 



^If one chooses 2mQ ~ 1, fo{z) — 5/4— ci and thus v^ax = 9/4, Vmin — 1/4. This gives cosh6' = 
(9/4-1/4) =5/4 and thus again agrees with the estimate on the basis of analyticity, 5/4 — ci = 0. This 
example shows that an ultra-local operator and an exponentially local operator are located very close to 
each other. 
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and also a large value for (2 — 2;)^'^+^. Thus our estimate of the solution is consistent for 
large k. 

We thus have 

1 TT TT 

1 — Ci~-(cos— ±isin— ). (5.21) 

If one sets Ci = cos(^o + i^) with a positive 5 and real we have 

Ci = cos(^o + i^) 

— cos ^0 cosh 5 — i sin Oq sinh 5 

1 i . TT , 

= 1 cos— =F-sin— ; (5.22) 

2 2A; + 12 2A; + 1 ^ ^ 

Since our equation (5.14) has real coefficients, we generally have a pair of solutions 

{60,6), {-00, S) (5.23) 

where 60 may be chosen to be positive; we thus explicitly solve only the positive solution 
in the following. For large k and small S, we have 

cos fc/n — 1 — - cos 



2 2A; + 1 ' 



sin^^o- Jl-{l-^cos^P 



2 ■ 



(5.24) 



From the relation 



we have 



1 TT 

sin 9q sinh 5 ~ - sin — (5.25) 

" 2 2A; + 1 ^ ^ 

1 . TT TT , 

sm — — ~J= • (5.26) 



2sin^o 2A; + 1 y^{2k 

From the above analysis of the Fourier transformation in (5.12), the crude estimate of 
the localization length, which is defined by 



is then given by 



\D{xi, 0)1 ~ ex.p[-x^6/a] = exp[-xi/L], (5.27) 



a/3 

L/a ~ — (2A; + 1). (5.28) 

71 



Generally, two (or more) solutions of ci = cos{9o+iS) with an equal imaginary part and 
different real parts super-impose the oscillating behavior on the exponential decay with a 
period of the order of the lattice spacing[19]. For the case k — (the overlap operator ), 
the closest singularity generally appears on the imaginary axis and no oscillating behavior 
occurs. 
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k 


Oo 


6 


L/a 


1 


1.69 


0.298 


3.35 


3 


1.62 


0.118 


8.47 


5 


1.60 


0.0737 


13.6 


7 


1.59 


0.0535 


18.7 


9 


1.59 


0.0421 


23.8 



Table 1: Localization length for k = 1,3, 5, 7, 9 

We also performed a numerical analysis of the values Oq + i6 in (5.22) and the local- 
ization length by solving the equation (5.14) numerically for k — 1,3, 5, 7, 9. (See Table 
1, where only the positive is shown). Applying the method of least square fit to these 
numerical data, we obtain the following relation 

L/a = 2.56A; + 0.8015. (5.29) 

This relation is not inconsistent with our crude analytical estimates (5.28). We reiterate 
that there exist two solutions with an equal imaginary part and different real parts and 
— , which explains the oscillating behavior super-imposed on the exponential decay in 
2-dimensional simulation[19]. 

Our estimates of the localization length (both of crude analytical estimate and numer- 
ical estimate) are much smaller than the general estimate in Ref.[9], which gives for the 
present one-dimensional example of H(2k+i), 

<=-h " = - 1 + 22-^<-«> (6.30) 

and thus a much larger localization length for A; ^ 1 

L/a ~ (2fc + 1)2^^ (5.31) 
In the notation of (3.44), our estimate (5.19) is based on 

/(cos 6)^1 + 2^''+\l - cos ef^+^ ~ (5.32) 

which suggests that 

Vma. = 2^(^'=+^) (5.33) 

for large k. For k ^ 1, our function f(cos9) in (5.32) has a very small slope near 
cos^ ~ 1, and thus the linear approximation gives a poor estimate of the solution of 
f{cos9) = 0. This particular property is not incorporated in the general Legendre expan- 
sion. In the general Legendre expansion, the huge maximum which appears in a specific 
sector of species doublers dominates the estimate of the localization length. We consider 
that this specific behavior of the integrand is the reason for the large discrepancy of our 
estimate from that of the general Legendre expansion [9]. The numerical simulation of 
the localization length in o? = 2 appears to support our estimate[19]. 
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6 Lower Bound for Interacting D^^^^^ D^^^^^ 

In this section we evaluate the lower bound for D^!'^^^^ D^^^^^^ and establish a finite 
locality domain of background gauge fields for the operator iJ^fc+i — H^2k+i) for a-ny finite 
k. This analysis is a direct extension of the analyses in Refs.[9][10]. 
We first define the directional parallel transporters [10]: 



= U^{x)ll}{x + /ia) = ^ Ui^(y)5a:+fia,y'(p(y). 



(6.1) 



Note that are unitary operators, and their norms are 1 due to the unitarity. Using 
and its Hermitian conjugate T^, D^r''^^^ is written as 



D 



(2fc+l) 
W 



n ^1 2A;+1 f„ N 2fe+l x2fe+l 

- ^{^7.(T,-rt)| +|^E(i-rt)(i-r,)| - (!^) .(6.2) 



We then have 



r,(2fc+l)tr,(2fc+l) 



rl ~|2(2fc+i) r-, ^ 

(' . ^ 2fc+l 

-'2[^i:i^-Ti){i-T,)j +1 



2(2fe+l) 



2k+l 



2fc+l 



(6.3) 



Here we set mo = 1, r = 1 and a = 1 . Noting that the third term in the above equation 
is a non-positive operator, we rewrite the above equation as follows, 



P,{2fc+l)t {2fe+l) 



2k+l 



2k+l 



2|E^(1-^/I)(1-^.)} 



2k+l 



+ 



+ 



r 1 1 

|i:7(i-r^)(i-rj(i-r;)(i-r„) 

Am 



2fc+l 
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- {e i(i - Tl){i - T,){1 - rt)(i - 
r 1 1 r 1 



2A;+1 



'1 

2fc+l' 



(6.4) 



Prom now on we refer to the expressions inside the four large square brackets in (6.4) as 
El, E2, Es, E4, respectively, and evaluate their norms in sequence. 

6.1 El 

We first define A = ((1 - T^,f + (1 - T^y) , B = E^(l - 7j)(l - T^)- Using A and 
B, El is written as 



El ^ [-T.\iT,-Tir^ +[T.li^-Tl){l-T,){l-Tl){l-T,)^ 

( -, ^ 2fe+l 

-2|E2(l-^i)(l-^^)| 

/-In 2fe+l /I ~' " 

^ {(_^ + 2S)2'=+^ + (A + 2S)2'=+^}-2f^ 



2fc+l 



2fe+l 

4; 



2( 2^^^^ + • • • + 2^^^^j + 2% 2A^''-^B^ + ■ ■ ■ + 2^M^^~^ ) 



+ . . . + 2^^-^( 2^^E^^-^ + ■ ■ ■ + 2B"'-'A') 

2k + lCl 



(6.5) 



In the last expression of (6.5), the notation such as 2fc+iCi shows the number of terms. 
We also note that only the even power of A and the odd power of B appear in the 
last expression of (6.5). Prom this fact, all operators appearing in the above equation 
can be rewritten as a sum of non-negative operators such as 2'+^74"^S^'+^A"^ and the 
terms containing a commutator S]. To show that 2^'^^ B^^'^^ are non- negative 
operators, we rewrite them as 

2/+2^m^2«+l^m = 2^+"^ A"' B^ (^{1 - T^)(l - T^)j B^A"" 
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by noting ^ A and = B. 

Now in each term in the last expression of (6.5) we want to bring all B to the center of 
the powers of A, but performing this manipulation for all the operators is quite tedious. 
To cope with this, we first note that the commutators, which appear when all B are 
brought to the center of the powers of A, are proportional to the field strength and that 
we only need the upper bound to those terms containing field strength when we finally 

evaluate the lower bound for D^!'~^^^ D^^~^^\ We therefore evaluate Ei approximately as 
the number of terms multiplied to the term which produces the largest number of terms 
containing the commutator. (By this way, we estimate the coefficient of the right-hand 
side of (6.7) below at a value larger than the true value. But this is a safe operation 
in the estimate of the upper bound.) The terms, which produce the largest number of 
terms containing the commutator, are the terms that include B at the right-most side as 
in A^k-2i j^2i+i ^ obtain {k — l){2l + 1) terms containing the commutator. 

As an example, we evaluate the norm of the term containing the (2/ + 1)— th power of 
B. Using the triangle inequality and < we obtain 



^22l+2^^2k-2l^2l+l ^ . . . ^ ^2l+1^2k-2l^ _ {^''^^^^^C^l+lA^-' B^'^^ A^' 



non— negative operators 



< 2''+\k+iC2i+i{k - l){2l + l)\\A\\''-''-'\\BnA, B]\\, l^O,l,---k- XB.7) 
Now we impose the constraint on the commutators of and as follows 

\\[T„n]\\ = \\[T„n]\\ = \\[Tl,n]\\ ^\\l-U,4<e, ^ u. (6.8) 
Then the norms of \\A\\, \\B\\ and are evaluated by using triangle inequality as 



||^||<32, ||5||<16, < 96e. 



(6.9) 



We thus have 



/-IN 2k+l k-l 

El - ( I ) E(2''^'2^+iC2;+iA'=-'S^'+iA^-') 



non— negative operators 



^22'+ViC2;+i(^-0(2^ + l)||^||''-''-iS||^'||[AS]|| 



2k+l k-l 




k-l 



< 3 Yl 2k+iC2i+i(k - l){2l + l)Uh. 

1=0 



(6.10) 



6.2 E2 

Here we define A = 
the terms in 



E2 = {A + Bf^+^ - A^^+^ 



Since all 
(6.11) 
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contain B and the norm of B is bounded by a value proportional to e, we can evaluate 
the upper bound to the norm of E2. As in the previous subsection, noting that ||A|| < 4 
and < 6e, we evaluate 



\E2\\ = \\{A + Bf''+^ -A^'^+^W < {\\A\\ + WBllf'^^ - 

2k+l 



2k+l 



(6.12) 



1=1 



6.3 Es 

We first rewrite E3 as follows, 



+ (1 - Tt)(l - T.)(l - Tt)(l - T,)) - ^ i (t,[TI, + Tt] + Tt[T„ + TJ]) 



2fc+l 



-{1:^(2-7^.-^5^} 



2fc+l 



A'^^{B + C) + --- + {B + C)A 



2k 



A^^-\B + C)^ + • • • + (S + CfA 



2 A2k-l 



H h < 



A(S + C)2'= + • • • + (S + C)2M 

^ V ' 

2fc+lC'l 



2^+1*^2 

) + {B + Cf^^\ 



(6.13) 



where 



B ^ E^{(l-rM)(l-Tt)(l-Tt)(l-Tj) 

+ (l-T;)(l-T.)(l-Tt)(l-r,)}, 
C ^ -Y.UT,[Tl,T^ + T^,] + Tl[T„n + T:^ 



(6.14) 



(6.15) 
(6.16) 



Here we note that A and B are non-negative and hermitian operators and that the norm 
of C is bounded by a value proportional to e. We split E3 into the terms containing C 
and the terms not containing C, and evaluate their norms in sequence. 
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6.3.1 Terms Not Containing C 

We further split those terms, which do not contain C, into the terms with an even power 
of A and the terms with an odd power of A. We first evaluate the terms with an even 
power of A. We then examine 

1=0 n 

Since ^'=-'52«+i^fc-« ^ non-negative operator, we want to write all the terms as a sum 
of those non-negative operators and the terms containing \A. B\. Here we adopt the same 
procedure as used in E\ above. The terms, which produce the largest number of terms 
containing the commutator, are the terms that contain B at the right-most side as in 
j^2k-2i^2i+i j^Yl the terms in (6.17) put together produce {k — l){2l + 1) terms which 
contain the commutator. We thus have 

||^2fe-a^2Hi + . . . + B2i+ij^2k~2i _ ^^^^C2i+iA^-^B^^+^A''-' II 

^ V ' 

non— negative operators 

< 2k+iC2i+i{k - m + l)Pf ^-2'-iS||2'||[A B]\\. (6.18) 

Noting that the norms ofA, 5 and are bounded as pll < 16, ||5|| < 48, ||[^,5]|| < 

192e, respectively, we obtain 

II ^(^2fe-2i^2Hl ^ . . . ^ ^2l+1^2k-2l^ _ 2k+lC2l+lA^-^B''^+^A''-^ \\ 
1=0 1=0 ^ " 



non— negative operators 
k 

< J2^k+iC2i+i{k - l){2l + l)162'=-2^-i482'l92e. (6.19) 

Next we examine the terms with an odd power of A 

k 

Y^^ ^^2k~2l+1^2l ^2lj^2k-2l+ly (g_20) 

We use the same procedure as above. Since J3'^2A:-2/+i^i ^ non-negative operator and 
the terms, which produce the largest number of terms with the commutator, are the ones 
that contain B at the right-most side, we obtain after counting all the terms in (6.20) 

II j^^j^2k-2l+l-Q2l ^ . . . ^ ^2/^2fe-2Hl) _ ^ 2ik+lCaS'>l2*^-2'+^S'|| 
1=1 1=1 

< j:2k+iC2il{2k -21 + l)\\Af'-'W'''\\[AB]\\ 
1=1 

k 

< J22k+iC2il{2k -21 + l)162'=-^'482'-^192e. (6.21) 

1=1 
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6.3.2 Terms Containing C 

Noting that the norm of C is bounded as ||C|| < 6e, one can confirm that the norm of the 
terms containing C is bounded from above by 

2fc+l m-1 

|2fc+l— m \ ' II D II ' II H '"1— ' 



E/^ II /i||2fc+l— m \ ' II Dll'll/^H'' 

m=l 1=0 
2fc+l m-1 



TO=1 /=0 

6.3.3 Final Result of Es 

Collecting all the caluculations in this subsection, we have 
/ 



1=0 



1=1 



V 



non— negative operators 



J 



2k+l 



< 2 E 2k+iCi l{2k + 1 - /)162'=-'48'e 
1=1 

2k+l m-1 

+ E E 2k+lCmmW^+^-"%S\Qtr-K 
m=l 1=0 



6.4 ^4 

We write i?4 as 

where 



= _,^2fc+1^2fc+l ^ ^^2fc+1^2fc+l 



B 



E(2-r,-Tj). 



(6.22) 



(6.23) 



(6.24) 



(6.25) 



The commutator [^2A:+i^ j^2fc+ij jg -^j-j^^gj^ ^s a sum of {2k + 1)^ terms containing the 
commutator We then have 

II _^^2fe+i^2fc+i ^ -^2/0+1^2^+1 II < {2k + lf\\Af^\\B\\''^\\[A,B]\\. (6.26) 

Noting that the norms of A, i? and [y4, 5] are bounded as ||y4|| < 4, ||i?|| < 8, ||[y4, _B]|| < 
6-\/2e, respectively^], we obtain 

11^411 < 6v^(2A; + l)2l6'=64^e. (6.27) 



•^Note that [^4, B\ can be rewritten as 

[AB] = [iT-^(T,-rt),i^(2-T,-T;)] 
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6.5 Lower Bound for Interacting 



r,(2/c+i)t r,(2fc+i) 



Finally, we combine all the caluculations in this section and obtain the lower bound for 
interacting D^^^^'^ D^^^^ as 

/I ^ 2fc+l 

> 111 +(7) Y.if''-\k+.C,i^iA'-^B'^^^A^-^) 

^ V ' 

non— negative operators 
J=0 Z=l 



non— negative operators 

2fc+l 



1=0 

l-E'211 



\l=0 1=1 



> 1- 



k-l 



1=0 
1=1 

2k+l 

+2 E 2/c+ia /(2A; + 1 - /)162^-'48^e 

2fc+l m-1 

+ E E 2fe+iC„mG16''+'-'"48^(66)'"-' 

m=l 1=0 

+6V2{2k + l)^16'=64'=e] . (6.28) 

By taking A; = 0, for example, in the above inequality, we have 

IIDvk^DvkII > l-6(2 + v^)e. (6.29) 

This inequality naturaly agrees with that derived by Neuberger[10] for the Wilson-Dirac 
operator, which corresponds to k = 0. The original result of Hernandez, Jansen and 
Liischer[9] gives a lower bound 1 — 30e. 

and that (7^ ± ± -fu)'' =2iov ji^v. 
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Taking k — l,for example, we obtain approximately 

> l-2xl0^e. (6.30) 

It is observed that the present estimate of the locality domain of gauge field strength for 
H{2k+i) gives a very small value. However, we have set all the non-negative operators in 

DyJ^ at in our estimate, as in (6.28), and those terms could be just as large as 
the coefficients of e, which are estimated by conservative huge upper bounds. We thus 
expect that the actual locahty domain of gauge field strength could be much larger. 

6.6 Upper Bound for D%^^^'^^ D%^^^^ 

We next evaluate the upper bound for D^j'^^^'^ D^l'^^\ Using the triangle inequality and 
||/li?|| < ||/1|| we can evaluate the upper bound for d'^^^^ as, 

\\Df^^\\ < i + 42'=+i + 8^"+\ (6.31) 
and then the upper bound for D^^'^^^^ D^^^^^ is evaluated as 

\\Df+^^^Df+^\\ < (l + 42^^+1 + 82'=+i)'. (6.32) 

6.7 Locality of Interacting H = a^y^D 

Our analysis of D^'''^^^^ D^'''^^^ estabhshes the locality domain of the operator 

i^(2.+i) - {7,aDr^' = + D^^'^'^ J ]. (6.33) 

for any finite k. Since this operator appears in the generalized chiral operator 

= 75 - i?(2fc+i) (6.34) 

our analysis provides a basis for the cohomological analyses of chiral gauge theory such as 
in [15] [16]. Note that defines the index TrT^ = n+ — n_. Our analysis in this section, 
which is regarded as an extension of chiral anomaly calculation, is naturally consistent 
with the analysis of anomaly in [13]. 

We have not established the locality domain of gauge fields for the operator 

H = a^,D = (//(2.+i))^/(^^+^) (6.35) 

explicitly in the interacting case. We however believe that H is local with approximately 
the same locality domain of gauge fields as for H(2k+i)- The reasons for this expectation 
are two-fold: The first is that we have established the locality of free H with a well-defined 
mass gap, which agrees with the mass gap for i/(2fc+i)- The second reason is a general- 
ization of (5.13). A multiplication of exponentially decaying operators also produces an 
exponentially decaying operator with the same localization length up to a possible poly- 
nomial prefactor. This suggests that the operation (6.35) also preserves the property of 
exponential decay for any finite k in the interacting case also. 
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7 Discussion and Conclusion 



We have discussed the locahty properties of the general class of lattice Dirac operators 
defined by (1.1), which satisfy the index theorem for any finite k. We first estabhshed the 
locality of all these operators at the free fermion level. We also presented a crude estimate 
of the localization length for these operators: The localization length approximately in- 
creases as 2k + 1, which is much shorter than the estimate on the basis of general Legendre 
expansion. It is clear that the operator, which is local at the free fermion level with a 
well-defined "mass gap" , is also local for sufficiently small and smooth background gauge 
fields. We in fact estabhshed a finite locahty domain of gauge fields for H{2k+i) for any fi- 
nite k by extending the analyses in [9] [10]. This locality is sufficient for the cohomological 
analyses of non-Abelian chiral gauge theory [15] [16], for example; note that our operators 
satisfy the exact index theorem on the lattice. We also argued that the locality domain 
of gauge fields for H itself is approximately the same as for i?(2fc+i) . From a view point 
of the regularization of continuum gauge theory, our operators with any finite k provide 
a satisfactory lattice regularization. However, our (conservative) estimate of the locality 
domain (6.28) is quite small. From a view point of numerical simulation, we thus need to 
perform a detailed numerical study of the locality domain and see if the actual locality 
domain is much larger, as we naively expect. Once a larger locality domain is established 
by numerical analyses, our general operators may be used for an explicit construction of 
chiral theories [15][16][20]-[24], for example. 

The present analysis of locality properties is also closely related to the perturbative 
evaluation of chiral U{1) anomaly for |agi74^| <^ 1 with a fixed lattice spacing, and the 
independence of anomaly coefficient for a small variation of the parameters r and mo [13]. 
See the similar analyses of anomaly for the overlap operator [24], which corresponds to 
A; = 0. 

As for the practical implications of our general operators, these operators exhibit prop- 
erties somewhat analogous to those of the "perfect action" of Hasenfratz et al. [26] ; both 
of these operators exhibit better chiral properties closer to those of the continuum oper- 
ator. Also the spectrum of the opeators with > is closer to that of the continuum 
operator in the sense that the small eigenvalues of D accumulate along the imaginary 
axis[17] compared to the standard overlap operator, for which the eigenvalues of D draw 
a circle in the complex eigenvalue plane. Contrary to an approximate solution to the 
perfect action[26], the explicit construction of our operators satisfies the index theorem 
exactly for all finite k. It may be interesting to apply our general operator with A; = 1, for 
example, to practical QCD simulations, since our operators preserve all the good chiral 
properties[27][28] of the overlap operator. 

One of us (K.F.) thanks T-W. Chiu for helpful correspodences about the numeri- 
cal analyses, which were essential for the present study, and Y. Matsuo for a clarifying 
comment on the Fourier transform of an analytic function. 
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